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PREFACE 



The purpose of this research was to find solutions of the diffe- 
rential equation of vorticity in gradient flew and to investigate the 
application of such solutions to the forecasting of upper air trajec- 
tories . 

This work was conducted at the U. S. Naval Postgraduate School, 
Monterey, California, during the period December 1950 to June 1951* 

I wish to express my appreciation to Professor W. D. Duthie for 
advice and guidance throughout. 
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I. INTRODUCTION 



The possibility of making aoourate mechanical forecasts of weather 
phenomena by stra ightf onward solutions of the basic equations governing 
atmospherio motions has long intrigued meteorologists. Theoretioally, 
such a solution to the problem is possible but there are so many varied 
and complex forces at work in the atmosphere that all attempts so far to 
effeot a complete solution to these equations have met with failure* 

One of the more important difficulties in the way of reaching satis- 
factory solutions has been the non-linearity of the equations and no 
standard methods are available to integrate nonlinear partial differential 
equations. Therefore the introduction of the so-called n perturbation" 
method was made which consists essentially in linearizing the equations by 
assuming meteorological quantities having basis values unchanging with time, 
plus small perturbation values whose second order terms oan be negleoted* 

One approach to quantitative forecasting has been the linearization 
and solution of the absolute vortioity equation. This was introduced by 
Rossby flOj , who obtained a solution for the motion of sinusoidal waves 
of infinite lateral extent in a horizontal plane and found that the velooity 
of propagation of such waves is given by the well known formula 



where C is the velocity of the waves toward the East , 1/ is average speed of 
the westerly current, L. is wave length and fl is the meridional rate of change 
of the Coriolis parameter. A considerable literature has grown out of the 
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interpretation of this formula on the classification of large scale 
atmospheric systems in terms of the so-called tonal index, •which is a 
measure of the strength of the zonal ourrent. 

Haurwitz extended Rossby' s method and obtained solutions for the 
motion of waves of finite lateral extent on a horizontal plane [5] and 
on a sphere [6j . Craig [ l] obtained solutions of the vortioity equa- 
tions in complete form without linearization for the plane and the sphere. 
These solutions differed from those of Eaurwitz by the absence of a zonal 
velocity term* Neamtan [ 7] then made a treatment of the vortioity in 
the same manner as Craig, viz by the use of stream functions and obtained 
identical solutions to those of Haurwitz. He showed that Craig's solutions 
were incomplete, causing absence of the zonal index term# 

Forsythe [ 3] developed formulas for the speed of propagation of waves 
with changing shape by use of the scalar relative vortioity as an identifiable 
property. However, he did not make any test on his formulas to determine 
their usefulness in practical forecasting. 

Ros sby and oo-workers [11] supplemented his initial formula based on 
the conservation of absolute vortioity with a more general teohnique which 
would be applicable to arbitrary initial streamline pattern so that the dif- 
ficulty of defining a prevailing wave length would be circumvented. 

All of the above solutions of the vortioity equations were predicated on 
the assumption of fr iotionless, non-divergent, and autobar otropic flow. 

i 

i ‘A * 

Further, in order to obtain a solution for relative vortioity which could be 
useful in forecasting trajectories it was necessary for Rossby Cll] to as- 
sume no horizontal shear of the wind and stationary pressure systems. 
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In this treatment, a general solution of the non-linear differential 
equation of the vertical component of absolute vorticity will be obtained 
•without as sximing non-divergent flow or baro tropic conditions. A solution 
is obtained for (1) a surface of constant height, and (2) a constant pres- 
sure surface. A discussion of mean-value constants for the wind shear term 
and for the effect of moving pressure systems, by means of which the vorticity 
equation could be integrated without disregarding these two terms, is included. 
A chapter is devoted to examining this solution with regard to its adaptation 
to forecasting of upper-air flow pattern by the trajectory method. In addition, 
special forms of the vorticity equation resulting from various assumptions are 
developed and discussed. 

Gradient flow is assumed throughout this investigation. The use of 
gradient flow permits a special handling of the divergence term. The wind 
of course is not always gradient even in the free atmosphere but is olocely 
approximated by the gradient wind at elevations greater than 1000 meters 
above the ground. However, under conditions of rapidly changing pressure 
gradient this close relationship between observed wind and gradient wind 
is greatly modified, and is due to the fact that the motion is not undor 
balanced forces. Under these conditions there is a velocity component 
along the isallobaric gradient. Consequently the use of gradient flow in 
the vorticity equation would generally be in error. 
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II. SOLUTION OF THE VORTICITI EQUATION FOR GRADIENT FLOW 



In obtaining solutions to the differential equation of the vertical 
component of absolute vortioity as few assumptions as possible will be 
employed. Initially, only the following assumptions are made* 

a. Friction is negligible 

b. Vertical velocities are negligible. 

Other restrictions will be added later to obtain solutions of particular 
oases. The general equation of vortioity to be solved, whioh is due to 
Bjerknes, is: 

y.\/+ 3 vr )/«, . (2 

M ~T J v 

It is noted that this equation consists of two terms, a term which 

/ ■» //jlM? j - 

represents the horizontal divergence of the gradient wind and a sole- 
noidal term which expresses the component of the geostrophio wind along 
the gradient of temperature* 



Since Va s V~P* &Z. , 



VJ\ ^(xr*Px£) 

( "t £ , J , J. I 



(2.2) 
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According to Taylor Tl2j the following expression for VyO)/ oan 



be derived: 



y./Oy'r - 'V-pxlk .y (?l + K'tV') , + C'h+ri s S) 



whioh, solved for V-V , yields 

V.y = - vpy/fe- 
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we obtain 



Now if we put this value of V- V in equation (2.2) and divido by ^ 

I . -*r 5T~ 

' 4V- ^ ^ ^ 

/ ^ 

y^t°< v^ ,e ‘ vCn+W)- (VPxK) _ (2<5J 



NOW V£>- Vf_ VT- 

/> p 



so that 



j__ £ _ 4. v/ £P_ 5 7r) I . ’ yptfc. \j (j\+kt v) •— "2J"! oi(VP< to) 

S«c*T~ P' r T Z 2 * 4 ) 



, iVf' ^ 

y, y p- o because raotors perpendicular* Henoe 

4- =_V.^X 4e <vP*/fe.. v(?>itivj- ¥T-~<6zf?L 

-5*. ~ ~r t* 



Now if we make the substitution 






/fe 



we obtain 

j-<±A = w. vrr 6+ kv ) ^ $+ ^ZS- CjuJ^A 
3, s* - - t A 



Rearranging and oolleoting terms 

)+P- 

"3C 



or 



rr(a t>^) Ji o> +kiv) v(t f «v)''j 

’_L_ ^ ^ W {[/l O* -tj+'d l vJL-r-i- vjLC'h+kiJ)} • 

c6f*- “ < ^ 



Integrating along the path of the partiole. 
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(2.5) 




wind, and the relation of curvature of streamlines to ourvature of trajeo- 



where the subsoript 0 indicates the initial values of the variable quantities. 

We note from equation (2.S) that on following a partiole along a path 
from an initial point to some other point on the path, the ratio of their 
absolute vortioities is equal to the ratio of their temperatures times the 
inverse ratio of terms involving curvature and shear times an exponential 
term. 

Equation (2.6) was derived without assuming either zero divergenoe or 
a barotropic atmosphere. Furthermore, no restrictions have yet been plaoed 
on the shear of the wind or ohanging streamlines* — 



tories 




Henoe equation (2*5) can be written, after 



integrating and rearranging terms. 







Now if we let 



and £ iz mean 
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If we assume a stationary pressure system, 9 as mos ^ 

investigations of vortioity, then equation (2 0 6) becomes 







( 2 . 8 ) 



1* Disoussion of Orders of Magnitude of the Mean Value Constant in 

Equation (2.8) for Stationary Systems. 

If we compare the relative order of magnitude of the terms contained in 
the exponential part of equation (2.8) some conclusions can be reached re- 
garding the value of the term^-J^ 80 that Q constant oan be assigned to 
this term which will represent its average value over most mid-latitude 
paths. With suoh a oonstant the exponential terms can then be integrated. 

An examination of 700 mb. oharts appears to demonstrate a radius of 
ourvature of most waves of the order of 600 miles. If we take an average 
value of 20 knots for the wind velocity the order of magnitude for the 
ourvature term will be 10"® sec"*. Estimating an order of magnitude for 
the shear term in a similar manner we use for an average shear a ohange of 
20 knots in a distance of 600 miles. This gives an order of magnitude of 
10"® seo"* for the shear term. An average value for the coriolis foroe 
term ^ would be the value at 45° latitude. At this latitude the ooriolis 
term is almost exactly equal to 10" 4 sec”*. .Therefore, the order of magni- 
tude for the ourvature and shear terms will closely approximate 10"® seo"* 
and for the ooriolis term will be 10"^ seo“*o It is seen from this quali- 
tative comparison that the ooriolis foroe is about ten times greater than 
the other two^terms in influencing ohanges in the absolute vortioity. 
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we obtain 



/ 9V 

TOien these values are placed in the expression > ^ 
the value of -'jj when the shear term is negative and the value of '/il 
when the shear is positive. 

Therefore the general solution (2.8) can now be written as 






A- tiA 

4 



( 2 .?) 



or 






where '/£ 

v 

The /_i- ) term will then be 



© 



‘ or (=^') * (2 differs 

from 2.0 by about yfo% Now, since the ratio of will usually be of 

the ordej' of unity the error involved in using simply instead of 

m will only be about 1%, 

Hence, by neglecting the exponential tern represented by tJ , equation 
(2.6) can be written simply as 

A , m '(Jd&k 






( 2 . 10 ) 



2. Discussion of Constants in the Equation for Moving Pressure Systems. 

All treatments of the vorticity equation appearing in the literature 
so far have ignored the effects of difference between the curvatures of 

streamlines and trajectories , and have assumed a steady state con- 

<c>t 

dition. This implies stationary pressure systems which of course is not 
the prevailing case in nature. The difficulty of handling mathe- 

matically has been responsible for neglecting this term. 
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Perhaps some oonolusions oan be drawn, however, from an examination 
of the range of variation of vO t-ki) so that a mean value oonstant 
oan be assigned to € with the result that equation (2*7) oan be inte- 
grated. 

From the relationship of trajeotory.ourvature and streamline curvature 
aooording to Petterssen C 8 , p. 225j ) ^ t***) , 

we can solve for v(ki-ki) so that 

= CC.S+ . 



Consider now a wave shaped system whioh is more or less symmetrical 
about a latitude oirole and moving eastward with speed C. Along the trough 
whioh is to the right of the path of the system Cos and is unity at 

bottom of trough. Along the ridge where is anticyolonio to left of 
path Cos If we ohoose average values of iiO.J for Cos ^ , the sign 

depending on whether along a trough or ridge, and take average values of 
1000 km and 20 kts for ^ and C respectively then s /(k* / 0 order 

of magnitude. 

From previous considerations of orders of magnitude for the shear and 
absolute vorticity we arrived at values of 10"^ seo"^ and lO*”^ seo"^ 
respectively for two terms. Therefore the final order of magnitude for 
will be two possible values j 




If the shear is negative, then £r^ — o-2 , and if shear is positive, 

6 . 
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Finally, substituting these two possible values in equation (2*7) 



we obtain 



and 




Therefore, under the conditions stated above the existence of positive 

r\ » 

or cyclonic wind shear cancels out the effect of the term and with 

anticyclonic or negative shear the correction term 6 i6 significant* The 
opposite effects would occur in the case of retrograde waves or in those 
parts of closed pressure systems where the wind blows in a direction 
opposite to the direction of motion of the system* 
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III. SPECIAL CASES OF THE GENERAL SOLUTION OF THE 
VORTICITY EQUATION FOR STATIONARY SYSTEMS 

Certain simplifying assumptions oan now be made with regard to the 
various terms of the general solution (2.8). The special oases whioh will 
be disoussed here aret 

1. Geostrophio Flow 

2. Wind Flow Parallel to Isotherms 

3* No Horizontal Wind Shear 

4. Combinations of the Above Three Cases. 

1. Case of Geostrophio Flow. 

Geostrophio flow prevails where there exist straight isobars or contours, 
henoe the ourvature term becomes zero and equation (2.8) beoomesj 



However, as demonstrated in Chapter II, the exponential term can be 
negleoted with only a small error involved so that for the geostrophio flow 
case equation (3.1) can be written as 



Equation (3.2) can be written as a quadratic and solved by the quadratic 
formula so that 




O 



(3.1) 









(3.2) 




(3.3) 



where 
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Equation ( 3*3) can now be expanded binomiallyj 




All terms beyond the first in this binomial expansion are of the order 



than for the shear term, oan be negleoted so that for the case of geostrophio 



Inasmuch as the absolute vortioity cannot be zero exoept for the very 
improbable oase of geostrophio flow with no shear and at the equator, the 
ohoioe of sign on the second term in (3»3) must be positive* 

For practical application on the weather ohart in determing how such 
flow would be affeoted by various types of movement, the relative vorticity 
is examined* For geostrophio flow, (3*2) can be written 



A particle moving northward would therefore undergo a deorease in 
relative vorticity or be turned anticyolonioally even though there may be 
a slight positive contribution from the temperature term* Southward flow 
would produoe the opposite effeot* 

A particle moving along a parallel of latitude would be influenced only 
by the change in temperature. If moving toward higher temperature, there 
will be an increase in relative vortioity and particle will ourve to the 



of 10"** or smaller and which, beoause of a muoh smaller order of magnitude 



flow the absolute vortioity can 



be expressed simply by = 




left. 
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2* Case of Wind Flow Parallel to Isotherms. 

The wind is blowing parallel to the isotherms in this case, which 
implies a barotropio atmosphere. Suoh conditions prevail when systems 
are thermally symmetrical* In this case equation (2*8) becomes 

- ( k? i V +'/)') o 

(kiV-f/)) (3.4) 



and the exponential term beoomes unity* Now, on following a partiole from 
initial point or path to some other point the absolute vortioity changes 
as the inverse ratio of their respective (fciV 4h) terms • 

If we now expand equation (3.4) a quadratic is obtained which can be 
written as X- X^ ~C , where C=(jf s i/+Jl (faV+T))* * the 



initial point values. Solving this equation by means of the quadratic 
formula gives 



(3*5) 



Equation (3*5) can now be expanded in a binomial expansion giving in a 
manner similar to Case It 




(3.6) 



The choioe of the positive sign for the second term is discussed under 
Case 1. 



Equation (3*4) can be now written as relative vortioity 

S= (ti V-t- a), + -s 

(ki /-/-/)) 
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3* Case of No Horizontal Wind Shear 



If the horirontal shear of the wind is either assumed to be tero or 
is negleoted then the general equation (2*8) simplifies to 



This oan be written as 





(3.8) 



We see then that for the oase of no shear the ratio of the absolute 
vorticity at two points is equal to the ratio of their respective absolute 
temperatures. It is also to be noted that this relationship whioh involves 
no assumptions as to divergence of the wind or barotropy differs from Rossby* s 
formula as applied to the forecasting of partiole trajectories only by the 
presence of the temperature term. 

If we express the vorticity as ^ and ohoose an inflection latitude, 

as did Rossby and oo-workers [ 11 ] , at which the ourvature is zero, then 
equation (3.8) oan be expressed as 



To 



(3.9) 



which corresponds to Rossby* s formula 

kiv — 7 )® - 7 > 

Equation (3.9) oould be applied in a manner entirely analogous to the 
technique used by Rossby in forecasting air partiole trajectories. The 
ooriolis force at any point from the inflection latitude oan be expressed 
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as a function of distance from such latitude so 



that J) ~ ft* + where ft 



is the rate of change of ooriolis foroe with latitude and is considered 
constant over the latitude range concerned* Then (3.9) can be written as 



If we express the curvature in terms of a second order differential 
equation of its y coordinate or distance from the inflection latitude we 
obtain after substituting this in (3.10) 



This differential equation oould perhaps be solved by elliptic integrals 
if a proper oonstant is chosen for the temperature term so that a value of \j 
coordinate would be obtained as a function of the intersection angle of 
particle path with inflection latitude* 

However, a qualitative interpretation of equation (3.10) oan be made 
showing effects on the relative vorticity with varying paths without actually 
solving the differential equation. If the particle is moving northward init- 
ially, the effect of coriolis force is to decrease the relative vorticity 
and increasing temperature tends also to increase the vorticity. However, 
now assuming decrease in temperature northward the temperature term in 
general is of much smaller magnitude than /3cj , henoe the particle curves 
anticyclonioally but to a slightly greater extent than the path resulting 
from coriolis considerations only. If the temperature increases northward 
along the path the effect is to decrease the effect of increasing coriolis 
foroe so that relative vorticity is decreased to a lesser extent than when 




(3.10) 
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the coriolis effect is considered alone. Therefore the path would curve 
less antioyclonioally than otherwise. 

For a southward moving ourrent the converse of all the above state- 
ments would hold* 

4* Combinations of the Three Preceding Cases, 

a. Geostrophio Flow with no Wind Shear. 

Under suoh conditions of wind flow equation (2.8) becomes simply 

3>. X' ? 



which can be simplified to 

_ 71 _ * 2 " 

7l* 

or the change of vortioity which is now due simply to change of 

ooriolis force )is directly proportional to the ohange of respective 

absolute temperatures. It is also interesting to note that the flow 

in this case implies a gradient of temperature from South to North, 

which of course, is not true in nature except under looal conditions. 

b. Geostrophic Flow with a Barotropio Atmosphere. 

In this case the general equation (2*8) can be written 

jf* 71 o 

3* 0 = ? 




0> *£)■ * 



c. A Barotropio Atmosphere with no TJind Shear. 
Under these conditions equation (2.8) becomes 

A- . OW* u- 

This can be written as 

(tSsY-f />) » . 
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This result is identioal in form to that used by Rossby and co- 
workers [ll~l in forecasting constant vortioity trajectories with the 
significant difference that no restrictions have been placed on divergence 
or convergence in this development. Therefore, the same technique as used 
by Rossby in forecasting partiole trajectories could be applied to the 
gradient wind using identioal relationships asj 

where ft - rate of change of coriolis force and (y) equals distance from 
an inflection latitude at which relative vorticity and curvature are zero. 
Solutions of the seoond order differential equation for the maximum y co- 
ordinate as a funotion of the intersection angle of the path with the in- 
flection latitude would be identioal to those of Rossby fill or Fultz C4]. 

The verification of forecasts of such computed trajectories should be 
more aoc urate than that of Rossby because of the added refinement of no 
restriction on divergence and convergence. 

d. Combination of all Three Speoial Cases — Geostrophio Flow, No 
Wind Shear, and Barotropy. 

These conditions would give a form for the vorticity equation thus 

Jl. 71 . 

•= 5 . "7 ' 

whioh gives the rather trivial result 

~ /lo 

However, absolute vorticity is again shown to be constant although only 
a function of latitude. Furthermore, zonal flow is implied by this 
equality of ooriolis terms. 
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IV. SOLUTION OF THE VORTICITY EQUATION FOR A CONSTANT 
PRESSURE SURFACE 



The equation of absolute vorticity along a surface of constant 
pressure is somewhat more simplified than along a oonstant height sur- 
face beoause the solenoid term is absent. This is not surprising since 
an isobar ic surfaoe obviously oannot be intersected by solenoids. 

Thus, equation (2.2) when applied to an isobario surface involves 
only a divergence term; 



TTe can now express this divergence in terms of ooriolis force, cur- 
vature and velocity in a manner entirely analogous to that for a constant 
height surface* 

Starting with the expression for the gradient wind; 







ait 



(4.1) 



y« 

" J 7\ 



which we oan express for an isobario surface as 




(4.2) 



Now multiplying by and rearranging; 



- f = y ^ ttvy 



solving for V gives 




Taking the divergence of both sides gives 



V. V z Vf * ■%> (h -fkiv) - (h+Hvj'v r - 3 Vr 
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whioh reduces to 



7 f .\/ =. V r *><&-V r (7H-tiv) 



-I 



(4.3) 



sinoe V r -V r <K-° ' 



It is seen that equation (4.3) gives the divergence of the gradient 
wind as the scalar produot of a vector parallel to the iso^ypses and the 
gradient of the quantity k*V) . 



Upon substituting this expression for the divergence in (4.1) and 

4 5*. w*l ^ M'”' 

rearranging we get i ^ • 

vM+mV. 

s. 3? y r 



Since can write, after performing 

indicated differential operations: 

-4- jy. y,(?+*iv) 

(7H-W) 

which can be expressed as the integral 



juOLf.)* T C \/.VjL(* +*’,*) . 



(4.4) 



Now since _ ]/ and Va." cLa-- 

d* 

we oan write (4.4) as 

^ aL ^*) - -j* 

whioh is readily integrable to 

L £- = f A , 

(Ti+tivX 

_ (7i + k?4v\ 

^ f **v~~ 



and 



(4.6) 
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Equation (4.5) is identical with (3.4) which was obtained for a 
constant level surface by assuming a barotropio atmosphere. However, in 
developing equation (4.5) no such assumption was necessary. For equation 
(4.5) to be applicable, however, the flow must be assumed to take place 
along an isobar io surfaoe. Therefore its application would be highly 
limited because of the improbability of isobar io flow. 

Now equation (4.5) can be rearranged as a quadratic and expanded 
binomially in a manner similar to equation (3.4) and whioh gives identical 



As a suggested extension to this work solutions of the vorticity 
equation oould be carried out in an entirely analogous manner for isentropic 
surfaoes and would have wider application than the solution for a constant 
pressure surface because of the greater prevalence of flow along isentropic 
surfaoes. 



results 
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